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ABSTRACT 



Coastal erosion, as inferred by measuring bluff recession is correlated with local 
wave height at twelve sites along the Monterey Bay coastline. Bluff recession rates are 
established by applying precise photo gramme trie techniques to a 44-year time series of 
aerial photographs. Wave heights are determined from the U'.S. Army Corps of 
Engineers Wave Information Studies spectral wave climatology, where deep water 
gravity waves are hindcast from historic wind fields at three-hour intervals from 
January 1956 to December 1975. The deep water spectra are refracted to shallow 
water spectra at a nominal depth of 4 m. An erosion model is developed for Monterey 
Bay where the erosion process is modelled as a non-linear function of the 4-m 
significant wave height 

R = [AHs 2 + BHs + C(Tide + 1.021 Is -ClilTtoe)]/Beach Slope 

The coefficients A, B, C are computed from a least squares regression of the modelled 
and observed recession rate values. The erosion model provides a reasonable 
representation of the erosion process in Monterey Bay, where the mean standard error 
between observed and modeled erosion rates is ±0.17 m/yr. Adjustment of the wave 
energy coefficient, A, allows tuning of the model for high and low wave energy 
locations. 
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I. EROSION PROCESSES IN MONTEREY BAY 



A. INTRODUCTION 

Each year, the energy generated in vigorous Pacific storms is released upon the 
Monterey Bay coastline. As the massive, long period waves break in the surf zone, 
their power is transformed into the forces that drive littoral transport and coastal 
erosion. These two processes scour and transport hundreds of thousands of tons of 
sand annually within the Bay. The popularity of Monterey Bay has brought pressure 
to develop its coastline; yet, past developments in this energetic zone, without proper 
regard for the coastal dynamics, have resulted in grave economic consequences. 



B. BACKGROUND 

Much of Monterey Bay area, particularly south of the Pajaro River, consists of 
low lying coastal plains bordered by foothills on the east, and active beaches on the 
west. This geomorphology is graphically depicted for South Santa Cruz County and 
North Monterey County in Figures 1.1, 1.2, and 1.3. Centuries of sedimentary 
discharge from the San Lorenzo, Pajaro and Salinas Rivers have provided the 
abundance of sand required to build the coastal plains and beaches. This discharge, 
combined with a prevailing northwesterly breeze, and generally southerly longshore 
transport, has generated expansive sand dunes along the inshore side of the beach. In 
recent geologic time, however, the trend is towards erosion and not accretion. Erosion 
occurs when powerful storm waves, coupled with high tides, are able to rush up the 
foreshore and scour away at the base of the sand dunes. The undermining of the cliff 
causes the face of the dunes to collapse and fall onto the berm. Successive waves 
entrain the slumped sediment and carry it out into the littoral transport zone where it 
is either deposited into a less energetic region, lost into the Monterey Canyon, 726,800 
m 3 /yr, Oradiwe (1986), or mined by one of Monterey's three sand companies, which 
collectively mine on the order of 413,100 m 3 /yr. Sand transport in Monterey Bay is 
shown in Figure 1.4. 
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Figure 1.1 South Santa Cruz County (Griggs, 1985). 
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Seawall 




Figure 1.2 North Monterey County (Griggs, 1985). 
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Figure 1.3 Central Monterey County (Griggs, 1985). 
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Figure 1.4 Littoral transport in Monterey Bay (Oradiwe, 1986). 
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C. THE NEED FOR COASTAL EROSION DATA 



The delightful climate and picturesque setting of Monterey and Santa Cruz 
Counties attract an ever expanding population. One manifestation of this population 
increase is mounting pressure to develop the coastal areas. Development of several 
areas along the Monterey Bay coastline without adequate knowledge of the coastal 
processes, has resulted in considerable monitary losses. 

Seacliff State Beach, at the northern extreme of our research area, is a well 
documented example of a location that experiences storm damage with cyclic 
regularity. The seawall at Seacliff State Beach was destroyed and rebuilt three times 
between 1940 and 1978. In 1980 a new seawall was constructed at a cost of 1.7 million 
dollars. The structure was engineered to last twenty years, but unfortunately, in the 
winter of 1983, storm waves and their associated flotsam battered the structure until 
one third of the seawall was completely destroyed, and heavy damage was inflicted on 
many nearby camping facilities. The estimated reconstruction costs were S740,000. 
(Griggs, 1985) 

Further south, the exclusive Pajaro Dunes condominium development was 
constructed in a region of active sand dunes. The development's location is one that 
experiences large scale episodic erosion followed by gradual accretion. The erosion 
occurs when the tidal stage is in phase with southwesterly storm waves. Three of these 
episodic erosion events occured between 1969 and 1983. The erosion was severe 
enough to cut the foredune back to the foundations of the development. Thousands of 
tons of riprap were placed in front of the development as a temporary measure during 
the 1983 storm season. The following summer, a permanent revetment fortress was 
built at a cost of several million dollars. (Griggs, 1985) 

The Moss Landing Marine Laboratory was constructed on a site that has some 
interannual variability; however, accretion is the dominant morphologic process 
occuring on longer time scales. The foundation for the laboratory was laid in 
mid-December, 1983; yet less than one month later, storm waves cut back the foredune 
and encroached upon the construction site. Again, tons of riprap were implanted as a 
temporary measure, and a permanent seawall was installed the subsequent summer. 

Perhaps the most dramatic and well documented case of episodic erosion has 
occurred at Stillwell Hall on the Fort Ord Military Reservation. Studies by Lima and 
Sklavidis (1984) found the average cliff recession rate to be 1.83 m/yr. _ This dramatic 
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erosion is consistent with the refractive focusing of wave energy to the area (Chapter 
III). The U.S. Army has recently spent S400,000 dollars as a temporary measure, while 
it contemplates whether to spend millions for permanent protection, move this soldier's 
recreation hall to a different site , or build a new one elsewhere. 

Lima and Sklavidis (1984) used precise, non-automated photogrammetric 
techniques to establish erosion rates at four locations in southern Monterey Bay: NTS 
Beach Lab, Monterey Sand Company, Phillips Petroleum property, and Stillwell Hall. 
Their results for the Phillips Petroleum property strongly corroborated the results of 
Thompson (1981), and their method of inferring coastal erosion through bluff recession 
(Thornton et. al., 1985) minimized the variance in photogrammetric measurements. 

Thompson (1981) viewed aerial photographs with a mirror steroscope to observe 
cliff toe recession at the Phillips Petroleum property in Monterey from 1939 to 1978. 
He concluded that the average recession rate was 0.6 m/yr. Moffitt (1968) used aerial 
photography to establish erosion rates in the vicinity of the Monterey Sand Company. 
Unfortunately, Moffitt referenced his data from the water line, the extreme variability 
of which resulted in rather large variance values. Nevertheless, he documented an 
average erosion rate of 2.1 m/yr. 

Lima and Sklavidis (1984) found cliff recession rates of 1.5 to 2.7 m/yr along a 
mile stretch of Sand City beach on either side of the Monterey Sand Company. They 
attribute this locally high erosion rate to the sand mining operations. 

The Monterey Beach Inn in Seaside and the Ocean House apartment complex in 
Monterey have both recently been threatened by major erosion events. The Monterey 
Beach Inn, previously named the Monterey Holiday Inn, constructed a permanent 
seawall to protect its investment, which at least thus far has been adequate. However, 
the Ocean House apartment management dumped five thousand tons of granite 
boulders onto the City beachfront adjacent to their complex as an emergency safegard 
measure. They are now required to remove the revetment from the public beach. 

The recent erosion events in Southern Monterey Bay may be natural phenomena 
in the annual and interannual variability of the coastline; or, as has been suggested by 
Griggs (1985), the increased erosion may be a result of sand mining in the Southern 
Bay. The reduction of sand availability caused by mining operations may inhibit the 
ability of large sand bars to form offshore during the winter. The annual formation of 
the offshore bars acts as a natural buffer against the powerful winter storm waves. The 
big waves break on these offshore sand bars and a majority of their energy is dissipated 



15 



within the surf zone. Without the bars, or if they are reduced in size, energetic winter 
waves can expend their force directly on the beachface, and provide the energy required 
for a major erosion event. 

The fact that several of the developed areas along Monterey Bay's coastline have 
suffered the severe consequences of building in the dynamic coastal zone should be 
taken as a forewarning. The increasing population of the Bay Area will bring more 
pressure to further develop this delicate area. It is clear that quantifiable erosion rate 
data and wave climatology must play a key role in site selection and structure 
engineering in order to prevent the inevitable destruction of improperly situated, or' 
engineered sites. 



D. OBJECTIVES 

The objective of this thesis is to develop a predictive model for coastal erosion in 
Monterey Bay. Erosion rates are inferred from bluff recession, which is measured using 
precise digital photogrammetric techniques from a 44-yr time series of aerial 
photographs. Reference points that define the beach profile must be chosen to insure 
both sharp definition of the reference on an aerial photograph, and accurate depiction 
of permanent, long-term erosion. The variability of the foreshore profile caused by 
seasonal sediment transport renders it unsuitable as a measurement reference. 
Similarily, cliff material that has fallen onto the cliff toe makes that potential reference 
difficult to define, and generally unreliable. The sharp stereo presentation of the cliff 
top however, offers an ideal profile, and one that is not subject to short time-scale 
variability. Consequently, changes in the cliff top profile reflect permanent recession. 
Twelve sites have been chosen from Monterey, north to La Selva Beach and bluff 
recession is measured along 400-to 1300-m coastal intervals centered on each of the 
twelve sites. The resultant erosion rate data is used along with a 20-yr spectral wave 
climatology and tidal information in the development of a coastal erosion model. 

A majority of the previous erosion studies that led to quantifiable results used 
imprecise methodology in photographic measurement. These inadequate procedures 
produced questionable results and unquantifiable errors. The precise photogrammetric 
technique used by Lima and Sklavidis (1984) produced results with acceptable 
confidence limits. The manual mode of the sfereocomparator that they used, however, 
induced operator measurement errors (reading the vernier), and was extremely time 
consuming. Hence only four sites were completed. 
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Our research employs digital equipment and techniques to minimize and make it 
easier to quantify operator and position measurement error. The variance in bluff 
recession data is minimized through numerical integration of the chosen coastal 
segments, and a comprehensive survey of the twelve coastal sites insures that the model 
will be representative of the erosion processes in Monterey Bay. 
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II. PHOTOGRAMMETRY 



A. BACKGROUND 

Precise photogrammetry is the currently recognized and practiced cartographic 
method of producing accurate, cost effective, topographic maps and shoreline 
manuscripts for nautical charts. Its accuracy for measuring events occuring on long 
time scales, such as cliff recession, is far superior to the previously employed methods 
of inferring coastal erosion through comparison of historical maps and charts, or 
subjective "hand" measurements on aerial photographs. 

This research uses planimetric and precise metric photogrammetric techniques to 
establish cliff recession rates. Planimetric photogrammetry is employed to qualitatively 
survey the selected measurement sites. At this initial stage, the photogrammetrist 
views the photo model stereoscopically to select the scaling points, identify local photo 
control points, determine the clarity of the cliff line, and establish the overall 
servicability of the photographs. Once these points have been identified, a field survey 
is performed to obtain the ground distance between scaling points, the height o'f the 
berm, and the elevation of the cliff top (above MLLW). Accurate quantitative 
measurements are then made using a precise stereocomparator. 



B. SITE SELECTION 

In any erosion study, it is advantageous to establish the erosion rates for the 
entire coastal survey area; unfortunately, the use of such a method on large scale 
erosion problems becomes economically unfeasible. A viable alternative, however, is to 
choose a number of study sites in the survey area that are representative of the overall 
erosion processes. This thesis examines twelve sites in Monterey Bay from the City of 
Monterey north to La Selva beach. Erosion information has been calculated for eight 
of the twelve sites examined, and results from Lima and Sklavidis (1984) are used for 
the remaining four sites. The specific criteria for site selection are listed below: 



* Is the site of economic interest ? 

* Is the site of scientific interest ? 
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* Is the site of political interest ? 

* Does the site lend itself to the use of photogrammetric techniques ? 

* Will the site provide representative erosion rates ? 

The total number of sites had to be limited to what was felt could be accomplished 
within the given time constraints. Additionally, the sites had to be spaced along the 
Bay's coastline in such a manner as to provide representative coverage, and all sites 
had to be photogrammetrically usable. The geographic positions of each of the sites 
and their associated names are provided in table I. The actual study area consisted of 
a 500-to 1 300-m stretch of bluff in front of each site. The site locations are shown in 
Figure 2.1. 



C. POTENTIAL ERRORS ASSOCIATED WITH PHOTOGRAMMETRY 

The primary sources of error in photogrammetric data are caused by image 
displacement and scale variation. Specifically, photogrammetric errors can be classified 
into three basic catagories: relief displacement, tilt displacement, or terrain induced 
scale variations. 

Relief displacement in a vertical photograph occurs when an image is displaced 
from its actual position due to an elevation difference between the image and the 
photographic principal point. Images of points on photographs with significant 
topographic relief will be displaced radially outward and images at lower elevations 
radially inward from the photographic principal point, as shown in figure 2.2. 
Displacement has the effect of portraying distances between images at higher elevations 
to be greater than they actually are, and distances between images at lower elevations 
to be less than they actually are. Relief displacement can be quantified and is given by 

d = (r x h)/ H (2.1) 

where d is the radial displacement distance, r is the radial distance from the principal 
point to the image, h is the image elevation above datum, and H is the flying height. 

Tilt displacement occurs when an image is displaced from its actual position 
resulting from a tilt in the vertical camera angle at the time of photo exposure. On a 
vertical photograph, the nadir, n, and the principal point, o, are indeed the same point. 
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TABLE I 

SELECTED SURVEY SITES 


Station Name 


Latitude 


Longitude 


NPS Beach Lab 


36°36T0" N 


121 °52'25" W 


Phillips Petroleum 


36°36'36" N 


121 °51'50" W 


Sand City 


36°37'00" N 


121°51T5" W 


Monterey Sand Co. 


36°37'20" N 


121 °50'55" W 


South Fort Ord 


36°38'05" N 


121 °49'55" W 


Stillwell Hall 


36°39'40" N 


121 °49'17" W 


Marina 


36°41'58" N 


121 °48'32" W 


Rincon Beach 


36°46'30'' N 


121 °47'48" W 


Marine Lab 


36°47'35" N 


121 °47'32" W 


North Moss 


36°48'25" N 


121 °47'25" W 


Zmudowski Dunes 


36°49'40" N 


121 °47'50" W 


Monterey Bay Academy 


36°54'26" N 


121 °50'55" W 
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However, when a camera is tilted at the time of exposure, its bore sight will describe 
the photograph's principal point, and a plumb line from the camera's optical center 
will describe its nadir, obviously now, not the same point. This distinction is important 
because image displacement on a tilted photograph will occur along radials of the nadir 
point, and not the principal point as was the case with relief displacement. 
Photographs with tilt displacement can be rectified by projecting their coordinates onto 
a new plane that is parallel to the map plane, or datum. The coordinates projected 
onto the new plane then comprise the equivalent vertical photograph. Tilt information 
was available for all of the photographs selected for this research, and only 
photographs with less than three degrees of tilt were used to establish clilf recession 
rates. The image displacement resulting from a 3° tilt of the horizontal flight line axis 
is less than the photocomparator's resolution ability; thus, all photographs employed in 
this thesis were considered vertical. The geometry of tilt displacement is given in 
Figure 2.3. 

Terrain-induced scale variation is another effect of vertical relief on a 
two-dimensional photograph. A manifestation of scale variation is that higher 
elevations will appear on the photograph at larger ’scales and lower elevations will 
appear at smaller scales (assuming a constant flying height). This variation is shown in 
Figure 2.4. In a photogrammetric mission, scale variation can be decreased by 
increasing the flying height and camera focal length. 

There are two procedures for minimizing the effect of the aforementioned errors 
in a photogrammetric data base. The first procedure is to apply point specific tilt, 
relief, and scale correctors to each measurement taken on the photograph. 
Unfortunately, hundreds of measurements are typically taken on each photograph, and 
this process is slow and laborious. The second procedure, and the method adopted by 
this research, is to choose ground control, and scaling points that are both, as close as 
possible to the measurement region, and at nearly the same elevation. Proper choice of 
the control and scaling points insures that the photogrammetric measurements are 
made over distances close to the baseline, and that the photoscale is highly accurate for 
the measurement region. Photoscale is given by 



S = (ab)/(AB) 



( 2 . 2 ) 
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Figure 2.2 Relief Displacement (MofTitt, 1980). 

where AB is the measured ground -distance from the field survey, and ab is the 
photodistance measured with a photocomparator. A more detailed discussion of 
photogrammetric errors and their correctors can be found in MofTitt (1980). 



D. INSTRUMENT ACCURACY AND SET UP 

Photogrammetric measurements were made on a modified Zeiss stereocomparator 
interfaced to an IBM Personal Computer (IBM-PC). A schematic of the Zeiss 
instrument is provided in Figure 2.5. The digital modification was accomplished by 
placing Teledyne-Gurley DMC-400 encoders on the stereocomparator's X, Y, and 
X-parallax adjustment controls. Each encoder outputs 8000 counts per one revolution 
of the comparator's controlling screws. Horizontal coordinates can be read from the 
instrument's verniers to an accuracy of 0.01 mm, and the X - parallax control can be 
read to 0.001 mm (Hallert, 1960). For calibration, numerous readings w r ere taken 
along each of the instrument's three verniers, and paired with their associated encoder 
counts. A linear regression was then performed for each of the three paired data sets 
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Figure 2.3 Tilt Displacement (Moffitt, 1980). 

to provide encoder count to horizontal distance transfer functions. These functions 
were evaluated for uncertainty by first calculating the standard deviation, or standard 
error for each axis, which were found to be: <r Y = 0.009 mm, <7 = 0.024 mm. The 

a y 

RMS radius of the error circle, or position error, d r , for an uncorrelated, orthogonal 
coordinate system is given by Heinzen (1980) as 

<t r = (<T 2 x + c 2 y )'/2 (2.3) 

and is found to be 0.022 mm . This photo position error can be converted to ground 
position error by multiplying d f by the photoscale. For example, the ground position 
error on a 1:12,000 scale photograph would be two tenths of a meter. 

Stereo vision presents the sharpest possible presentation of the standard 9" x 9" 
photographic data. A pair of overlapping photographs are placed on plates 1 and 2 of 
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Figure 2.4 Scale Variation (Moflitt, 19S0). 

the photocomparator (Fig. 2.5), each of which may be rotated in a plane normal to the 
vertical axis of the photograph. The image on plate 1 is rotated until the X-fiducials 
are parallel to the flight line (X-axis of the photocomparator). The image on plate 2 is 
also rotated to align the X-fiducials with the flight path, but additionally, a Y-parallax 
adjustment must be made to bring the flight line on image 2 into coincidence with the 
flight line on image 1. The Y-parallax will have been completely removed once the two 
flight lines are in coincidence, and the operator will be able to see two marks in the 
stereo viewer. One mark, on image 1, will be stationary, the other mark will appear to 
float above the photograph ( Figure 2.6). The X-parallax or elevation, control is then 
adjusted until the floating mark becomes coincident with the fixed mark. The 
conjugate points on the two photagraphs are now viewed as exactly the same point 
and precise stereo measurements can be made. 



E. MEASUREMENT PROCEDURE 

Several points are selected that appear on each stero model in a particular 
photographic time series. These points are: a coordinate system origin, O; a baseline 
definition point, B; and numerical integration limits, Lj and L u (Figure 2.7). Annual 
and interannual variability of the sand dunes render it impossible to repeatedly locate 
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Figure 2.5 Zeiss Sterophotocomparator. 

the integration limits by sight on a photograph; instead, these points are accurately 
located by the intersection of the clifT edge and a line normal to the clifF that has been 
defined by two reliable points, I Pj and 1P 2 discernable on the photographs. The same 
scaling points are typically selected for each stereo model. This procedure is not 
always possible, however, and additional scaling points sometimes have to be chosen. 
Such additional points must be supported by supplimental field work. 

Digital coordinates are obtained through employment of an interactive program 
which is designed to communicate between the comparator's encoders and the 
IBM-PC. The program, PHOTODAT, accepts commands from the 1BM-PC keyboard 
and writes the desired coordinates into a station data file. PHOTODAT is listed in 
Appendix A. The digital data acquisition procedure is very specific and requires the 
origin of the photographic coordinate system be defined first. This step is 
accomplished by moving the comparator mark to the designated position, O, and 
entering it. The operator must then move the mark to the baseline definition point, B, 
and enter its position. The X and Y values of point B define the coordinate system 
rotation angle a given by 
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Figure 2.6 X-Parallax and the Floating Mark, 
o = arctan(Y/X) * (2.4) 

All points entered subsequent to the baseline definition point are projected onto a new 
coordinate system via the coordinate system rotation matrix 



x r 




cosa sina 




~x“ 


Y 1 


= 


-sina cosa_ 




_Y_ 



(2.5) 



The scaling parameters are the next points entered, and the first points to be projected 
onto the new coordinate system. Once the photographic scale has been accurately 
defined, blufF measurenents can begin. The clifT edge is outlined by a dense network of 
points entered monotonically along the clifT from point Lj to L u> The actual cliff 
outline is approximated by line segments whose lengths vary slightly with topography 
but are typically less than 0.25 mm at photoscale. A trapazoidal rule integration was 
performed to determine the area between the cliff edge and the baseline using 

B( Y ; + Y i+i)/2 X |X i+ i -Xjl] (2.6) 

t'l 
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Areas for adjacent years in the time series were subtracted to provide the actual eroded 
area, and that figure was multiplied by the mean cliff elevation to determine the eroded 
volume. The actual coordinate system rotation, scaling, and numerical integration are 
performed by an interactive program, WRKPROG (Appendix B). WRKPROG 
retrieves raw data from the PHOTODAT files, accepts inputs from field surveys, and 
executes the numerical computations. 



F. FIELD SURVEY 

A separate field survey was performed for each of the eight erosion sites. One 
task of the field party was to determine the authenticity of the photo control and 
scaling points. If there was any doubt as to whether a particular reference point had 
been moved or changed sometime during the 44-year photographic time series, that 
point was rejected, and another selected. Additionally, accurate distances had to be 
measured between the ground scaling points to determine the precise photographic 
scale. This operation was accomplished using a Keuffel and Esser, Rangcr-lV, laser 
distance measuring device. The K&E Ranger-IV has an accuracy of ±5 mm 4- 2 ppm 
(Keuffel and Esser, 1975). At typical measurement distances of 800 m, this uncertainty 
translates to ± 5 * 10' 4 mm on a 1:12,000 scale photograph. A standard deviation of 
±5 x 10' 4 mm is one full order of magnitude less than the Zeiss photocomparator's 
resolving power; consequently, no accuracy was sacrificed through use of the 
Ranger-IV. The scaling distances at each station were measured ten times and 
averaged to insure the best possible results. Temperature, pressure and humidity were 
observed at each measurement site, and the measured distance was corrected for 
atmospheric refraction effects. 

The elevation of the cliff toe was measured at 75-m intervals along each site's 
beachfront. The elevation measurement was observed using a surveyor's level and a 
Philadelphia rod. The actual cliff toe elevation was leveled from the mean water line, 
and the measurements were subsequently reduced to datum (MLLYV) using tidal 
correctors. 
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G. EROSION DATA PRESENTATION 



Erosion results are presented for eight sites along Monterey Bay: Sand City, 
South Fort Ord, Marina, Rincon Beach, Moss Landing Marine Lab, North Moss 
Landing, Zmudowski Dunes, and Monterey Bay Academy. 

Typically, 75 to 150 points are digitized from each stereo model. Integrated areas 
for adjacent photo models in the time series are subtracted, and the difference is the 
recessed area for the intervening time period. The recessed area is divided by the 
longshore observation distance (length of the digitized bluff), and that value is divided 
by the number of years between the photo models. The result is a cliff recession rate 
in meters per year. Negative values in the tables below indicate recession, whereas 
positive values indicate accretion. A standard deviation is calculated for the recession 
rates at each site, and that value is stated in the respective site's result table. The 
number of standard deviations by which a particular recession rate drifts away from the 
mean value is listed in the tables' "|xj - x|/cr " columnns, and provides a measure of 
temporal variability. Assuming episodic recession, a high value in the sigma column 
would signal an erosion event; however, it is important to also consider the length of 
the intervening time period. Long time intervals will tend to alias erosion events, and 
conversely, erosion events occuring within a short time period will be highlighted. 
Average erosion is inferred by multiplying the recession rate by the average cliff 
elevation. Not surprisingly, areas of refractive wave focusing resulting in energy 
convergence experienced the highest recession rates; whereas sites located at the head 
of Monterey Canyon, where there is refractive divergence, experienced long-term 
accretion. It is noteworthy that all sites experienced at least one significant erosion 
event between 1980 and 1984. 

Sand City is a site that does not have a well defined coastal bluff. The bluff was 
approximated by digitizing the crest of the seaward-most sand dune. While recession is 
the dominate process occurring over the 44-year time series, significant accretion 
occurred between 1940 and 1956. This point is interresting, because - whether this 
growth event was a manifestation of sand dune variability, or bonifide accretion - the 
process made an abrupt change late in the time period. In the years following 1966 
Sand City displayed recession rates comparable to those of Fort Ord, where high wave 
energy drives the highest recession rates recorded in Monterey Bay. The high recession 
rates observed at Sand City are not corroborated by the wave energy studies discussed 
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in Chapter III. The recession results at Sand City most likely reflect the presence of 
the sand mining industry which is most prevalent in that region of the Bay (Table II 
and Figure 2.8). 

The highest recession rates we observed were at South Fort Ord. The recession 
is depicted in Table III and Figure 2.9, where the average bluff receded fifty meters 
from 1940 to 1984. This high rate is consistent with the refractive wave energy 
focusing which appears to be a maximum in this region. 

Marina is also a site where high wave energy produces high recession rates. 
Additionally, numerous observations suggest that Marina is an area of high wind 
energy, perhaps due to geographic focusing of the seabreeze. The high wind energy 
and the variable nature of the non-solidified coastal bluff cause high variance in the 
observed recession rates, and may, in part, be responsible for the anomolous accretion 
observation in the 1970 time frame (Table IV and Figure 2.10). 

Station Rincon Beach at the northern "edge" of the South Bay's energy 
convergence zone displays continuous recession from 1940 to 1984. The inference of a 
convergence zone boundary is made as wave energy increases south of Rincon Beach, 
and stations just north of Rincon exhibit long-term accretion. The recession at Rincon 
Beach is seen in Table V and Figure 2.1 1. 

Accretion at 'station Marine Lab is evidenced by a shift in the entrance of 
Elkhorn Slough which occurred sometime between 1940 and 1966. This shift is clearly 
depicted in Figure 2.12. Figure 2.13 omits the slough entrance to present the 
recession/accretion events as they occurred after 1966. The numerical values are seen 
in Table VI. 

The low wave energy caused by refractive divergence at the head of Monterey 
Canyon is responsible for littoral deposition and the resultant accretion at station 
North Moss Landing (Table VII and Figure 2.14). 

Station Zmudowski Dunes appears to be the northern conjugate of Rincon 
Beach. It is a location where higher wave energy produces continuous recession, yet 
only 1.5 miles to the south, lower wave energy allows the accretion process to 
dominate at station North Moss Landing. Figure 2.15 graphically depicts the 
recession, and the recession rates are provided in Table VIII. The variability of 
McClusky Slough is evident in figure 2.15 and is omitted to present Figure 2.16. 

Recession is the dominate process occurring at Monterey Bay Academy between 
1956 and 1984. The wave energy focused to the north of Monterey Canyon is typically 
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not as high as energy which has been focused to the south. This point is almost 
intuitive, as a majority of the winter swells originate from extratropical cyclones in the 
Gulf of Alaska. Accordingly, the -recession rates at Monterey Bay Academy are less 
than those at its southern conjugate station South Fort Ord. The recession at 
Monterey Bay Academy is presented in Table IX and Figure 2.17. 
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TABLE II 

EROSION RATE RESULTS FOR SAND CITY 



STATION NAME : Sand City 
LATITUDE : 36°37'00" N 
LONGITUDE : 121°51'15" W 

MEAN CLIFF TOE ELEVATION (above MLLW) : 4.53 M. 
MEAN CLIFF TOP ELEVATION (above MLLW) : 10.04 M. 
LONGSHORE OBSERVATION DISTANCE : 540.00 M. 
BEACH SLOPE : 1:40.00 
SIGMA : 1.00 M. 



TIME INTERVAL 


RECESSION 

[M/YR] 


|x r x| 

(T 


AVG. EROSION 
KM 3 /YR)/M] 


1940 - 1956 


+ 0.47 


1.44 


+ 4.72 


1956 - 1970 


-2.30 


1.33 


-23.09 


1970 - 1976 


-1.22 


0.25 


-12.25 


1980 - 1984 


-1.34 


0.37 


-13.45 


1940 - 1984 


- 0.97 




-9.74 
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Figure 2.8 Sand City Erosion Results. 
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TABLE III 

EROSION RATE RESULTS FOR SOUTH FORT ORD 



STATION NAME : South Fort Ord 
LATITUDE : 36°38'05" N 
LONGITUDE : 121°49'55" W 

MEAN CLIFF TOE ELEVATION (above MLLW) : 6.19 M. 
MEAN CLIFF TOP ELEVATION (above MLLW) : 35.39 M. 
LONGSHORE OBSERVATION DISTANCE : 692.00 M. 
BEACH SLOPE : 1:36.36 
SIGMA : 1.47 M 



TIME INTERVAL 


RECESSION 

[M/YR] 


|x r x| 

<7 


AVG. EROSION 
[(M 3 /YR)/M] 


1940 - 1956 


-1.09 


0.02 


-38.58 


1956 - 1966 


-0.01 


0.75 


-0.35 


1966 - 1980 


-1.13 


0.01 


-39.99 


1980 - 1984 


-3.85 


1.86 


-136.85 


1940 - 1984 


-1.12 




-39.64 
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Figure 2.9 South Fort Ord Lrosion Results. 
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TABLE IV 

EROSION RATE RESULTS FOR MARINA 



STATION NAME : Marina 
LATITUDE : 36°4U58" N 
LONGITUDE : 121 °48'32" W 

MEAN CLIFF TOE ELEVATION (above MLLW) : 4.92 M. 
MEAN CLIFF TOP ELEVATION (above MLLW) : 24.31 M. 
LONGSHORE OBSERVATION DISTANCE : 1283.00 M. 
BEACH SLOPE : 1:40.00 
SIGMA : 1.67 M. 



TIME INTERVAL 


RECESSION 

[M/YR] 


|x r x| 

<7 


AVG. EROSION 
[(M 3 /YR)/M] 


1940 - 1956 


-0.99 


0.36 


-14.83 


1956 - 1970 


-1.12 


0.44 


- 27.23 


1970 - 1980 


+ 2.30 


1.60 


+ 55.91 


1980 - 1984 


-2.12 


1.05 


-51.54 


1940 - 1984 


-0.38 




-9.24 
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Figure 2.10 Marina Erosion Results. 
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TABLE V 

EROSION RATE RESULTS FOR RINCON BEACH 



STATION NAME : Rincon Beach 
LATITUDE : 36°46'30" N 
LONGITUDE : 121°47'48" W 

MEAN CLIFF TOE ELEVATION (above MLLW) : 5.33 M. 
MEAN CLIFF TOP ELEVATION (above MLLW) : 8.41 M. 
LONGSHORE OBSERVATION DISTANCE : 1353.00 M. 
BEACH SLOPE : 1:69.56 
SIGMA : 1.26 M. 



TIME INTERVAL 


RECESSION 

[M/YR] 


|x r x| 

<J 


AVG. EROSION 
[(M 3 /YR)/M] 


1940 - 1956 


-0.82 


0.02 


-6.89 


1956 - 1966 


-0.20 


0.51 


-1.68 


1966 - 1980 


-0.63 


0.17 


-5.29 


1980 - 1984 


- 3.27 


1.93 


-20.44 


1940 - 1984 


-0.84 




-7.06 
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Figure 2.11 Rincon Beach Erosion Results. 
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TABLE VI 

EROSION RATE RESULTS FOR MARINE LAB 



STATION NAME : Marine Lab 
LATITUDE : 36°47'35" N 
LONGITUDE : 121°47'32" W 

MEAN CLIFF TOE ELEVATION (above MLLW) : 5.15 M. 
MEAN CLIFF TOP ELEVATION (above MLLW) : 8.93 M. 
LONGSHORE OBSERVATION DISTANCE : 506.00 M. 
BEACH SLOPE : 1:66.66 
SIGMA : 1.63 M. 



TIME INTERVAL 


RECESSION 

[M/YR] 


|Xj-x| 

(7 


AVG. EROSION 

[(m 3 /yr)/mj 


1940 - 1966 


+ 1.19 


0.45 


+ 10.63 


1966 - 1970 


-2.04 


1.53 


-18.22 


1970 - 1976 


+ 1.30 


0.52 


+ 11.61 


1976 - 1984 


-1.35 


1.10 


-12.05 


1940 - 1984 


+ 0.45 




+ 4.02 
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Figure 2.12 Marine Lab Erosion Results (Unsmoothed). 
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Figure 2.13 Marine Lab Erosion Results (smoothed). 
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TABLE VIJ 

EROSION RATE RESULTS FOR NORTH MOSS LANDING 



STATION NAME : North Moss Landing 
LATITUDE : 36°48'25" N 
LONGITUDE : 121°47'25" W 

MEAN CLIFF TOE ELEVATION (above MLLW) : 5.29 M. 
MEAN CLIFF TOP ELEVATION (above MLLW) : 9.81 M. 
LONGSHORE OBSERVATION DISTANCE : 900.00 M. 
BEACH SLOPE : 1:66.66 
SIGMA : 2.08 M. 



TIME INTERVAL 


RECESSION 

[M/YR] 


|x r x| 

<7 


AVG. EROSION 
l(M 3 /YR)/M] 


1940 - 1970 


+ 0.09 


0.03 


+ 0.88 


1970 - 1976 


+ 2.64 


1.26 


+ 25.99 


1976 - 1980 


-0.12 


0.07 


-1.47 


1980 - 1984 


-3.21 


1.56 


-31.78 


1940 - 1984 


+ 0.03 




+ 0.29 
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Figure 2.14 North Moss Landing Erosion Results. 
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TABLE VIII 

EROSION RATE RESULTS FOR ZMUDOWSKI DUNES 



STATION NAME : Zmudowski Dunes 
LATITUDE : 36°49'40'' N 
LONGITUDE : 121°47'50" W 

MEAN CLIFF TOE ELEVATION (above MLLW) : 3.26 M. 
MEAN CLIFF TOP ELEVATION (above MLLW) : 15.56 M. 
LONGSHORE OBSERVATION DISTANCE : 1050.00 M. 
BEACH SLOPE : 1:69.55 
SIGMA : 0.28 M. 



TIME INTERVAL RECESSION 

[M/YR] 


|x r xt 

<T 


AVG. EROSION 
[(M 3 /YR)/M] 


1940 - 1956 


-0.16 


1.28 


-2.48 


1956 - 1970 


-0.84 


1.14 


-13.07 


1970 - 1984 


-0.61 


0.32 


- 9.73 


1940 - 1984 


-0.52 




-8.09 



46 



ZM-UDOWSK I DUNES 




LD 

CM 



LD 

CO 



in 

r\j 

ro 



LD 

r\ 

CM 



LO 

CM 

CM 



(S33J.3W) NO I SS3033 



Figure 2.15 Zmudowski Dunes Erosion Results (Unsmoothed). 
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Figure 2.16 Zmudowski Dunes Erosion Results (Smoothed). 
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TABLE IX 

EROSION RATE RESULTS FOR MONTEREY BAY ACADEMY 



STATION NAME : Monterey Bay Academy 
LATITUDE : 36°54'26" N 
LONGITUDE : 121°50'55" W 

MEAN CLIFF TOE ELEVATION (above MLLW) : 4.36 M. 
MEAN CLIFF TOP ELEVATION (above MLLW) : 21.74 M. 
LONGSHORE OBSERVATION DISTANCE : 1220.00 M. 
BEACH SLOPE : 1:76.79 
SIGMA : 0.64 M. 



TIME INTERVAL 


RECESSION 

[M/YR] 


|x r x| 

d 


AVG. EROSION 
[(M 3 /YR)/M] 


1940 - 1956 


+ 0.18 


0.79 


+ 3.68 


1956 - 1970 


-1.33 


1.45 


-30.07 


1970 - 1984 


-0.03 


0.04 


-0.60 


1940 - 1984 


-0.33 




-6.62 
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Figure 2.17 Monterey Bay Academy Erosion Results. 
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111. WAVES AND TIDES 



A. WAVE CLIMATOLOGY 

An accurate understanding of wave climatology is paramount in any study of 
beach erosion. Wind waves provide the energy to drive the nearshore dynamics 
responsible for changes in coastal morphology. The wave climatology for Monterey 
Bay was calculated using the Wave Information Studies (WIS) deep water variance 
density spectra compiled by the U.S. Army Corps of Engineers. The WIS data is 
broken down into three phases. This study will be concerned with phases 1 and II 
only. 

In phase I, wind fields are derived from the Gridded Northern Hemisphere 
Pressure Fields of the Navy Fleet Numerical Oceanography Center (FNOC), and the 
North American Historical Weather Map Series of the National Climatic Center 
(NCC). The data from FNOC and NCC represent a twenty year sampling from 1 
January 1956 through 31 December 1975. A first guess wind field in regions of low 
streamline curvature is computed from quasi-geostrophic theory. Baroclinic effects 
caused by vertical temperature gradients (thermal wind) are added to the wind 
components and provide the basis for the observed ageostrophic flow. Wind fields for 
synoptic systems exhibiting strong streamline curvature, as typified by an extratropical 
cyclone, are derived from gradient wind theory which accounts for centrifugal 
acceleration in the balance of forces. Additionally, an isallobaric component of the 
wind field is estimated and included in rapidly developing systems where the rate of 
change in pressure with time cannot be ignored. The derived free atmospheric wind 
fields are then reduced to a boundary layer reference level using friction velocity scaling 
(Resio, Vincent, and Corson, 1982). Phase I wind information is computed with a six 
hour time step and interpolated to a 220 KM open ocean grid. In Phase II Surface 
Marine Observations and Airways Surface Observations, both available from NCC, are 
blended into the wind fields computed in phase I, and the modified wind fields are 
interpolated to a 55-km coastal grid. 

Wind waves are hindcast using a discrete spectral wave model developed by 
Resio (1981). Two-dimensional variance density spectra are generated over sixteen 
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directional bands at 0.01-Hz frequency bands from 0.03 Hz to 0.22 Hz. Discrete 
spectral models, in general, are representative for waves in their late stages of 
development; however, they exhibit significant difficulty simulating proper growth in a 
fetch or duration limited domain, or for waves in their early stages of development 
(Dexter, 1974; Resio and Vincent, 1979). Barnet (1968) proposed a more complete 
specification of the classical A + BE(f) source term as a possible solution to this 
problem, where A represents the Phillips (1957) linear, initial wave growth mechanism, 
B , the Miles (1957) exponential, continued wave growth mechanism, and E(f) is the 
energy density spectrum. Although Barnet's term was three times too low for the 
JONSWAP spectrum, his initiative in incorporating a non-linear input aided in the 
subsequent parameterization of non-linear wave-wave interaction used in the Resio 
(1981) model. The non-linear source term employed by Resio (1981) is given by 

G nI = Da 3 g 2 f m 4 exp[l - (f m /f) 4 ](f/f m ) 3 (3.1) 

From the JONSWAP data, a = 0.076X, where X is the nondimensional fetch 
X= gx/U 2 , D is a dimensionless constant, f m is the frequency of maximum energy, and 
U is the wind velocity. The linear Phillips mechanism was found to transfer more 
energy into the low frequencies on the spectrum's leading edge than into the mid-range 
frequencies. It was consequently removed from the Resio (1981) model and replaced 
with the G n j term in the source equation. The total source term for the Resio model 
is then given by the non-linear wave-wave interaction parameterization and the Miles 
exponential growth mechanism 

Source “ G„, + B E(f) V(0 - 9 0 ) (3.2) 

In Equation 3.2, \|/ is a nondimensional directional spreading function where 0 Q is the 
mean wind direction, and B is given by 

B = s f{(u cos0)c _1 - 0.9] (3.3) 

In Equation 3.3, s is the density ratio between the air and water, and c the phase 
speed, or celerity. 
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Clearly, with the Phillips mechanism missing, another wave generation scheme 
must be included to provide initial wave growth. This initialization, is accomplished 
with a local parametric model. The domain of the parametric model is limited by a 
fixed high frequency cut off which delineates the boundary between the discrete 
spectral and parametric domains. The peak frequency is calculated for each time step 
from equations presented by Hasselmann et. al. (1976). The energy developed in the 
parametric domain is transferred to the discrete spectral domain until the peak 
frequency reaches the high frequency cut-off. The parametric spectrum then shifts to 
an f 5 equlibrium form, and all growth is continued in the discrete spectral domain. 
The Resio (1981) model provides wave growth and decay rates that are consistent with 
the observed JONS WAP data while maintaining spectral shapes that are highly 
representative of observed spectra (Risio, 1981). In phase I, wave information is 
calculated on a 220-km grid and the spectra are archived at 64 locations in the Pacific 
as depicted in Figure 3.1. in phase II, wave information is calculated using the 
modified wind fields at 55-km grid spacing and the spectra are archived at 53 coastal 
locations as depicted in Figure 3.2. Note that the wind wave model is a northern 
hemispheric model, and it does not allow wave energy from the southern hemisphere to 
propagate across the model's boundary. 

B. WAVE REFRACTION 

Wave rays that approach non-parallel isobaths at an oblique angle will tend to 
bend toward the orthogonals of the isobaths as a result of wave refraction. 
Consequently, an area where rays come together is a region of energy convergence, and 
an area where rays separate, a region of energy divergence. In Monterey Bay the Moss 
Landing area at the head of Monterey Canyon provides a good example of energy 
divergence; conversely, Fort Ord typifies a region of ray focusing, or energy 
convergence. 

The effects of shoaling (changing group velocities) and refraction (convergence or 
divergence of variance density) are applied to the WIS deep water spectra to provide 
shallow water and surf zone information. The variance density transformation is 
performed using spectral ray theory. In spectral ray theory, one assumes that energy, 
or variance, that is associated with a particular frequency and directional band will 
remain in that band throughout the transformation, i.e. constant energy flux 
(Thornton, 1983). Hence the shallow water variance density spectra S^(f,0) are related 
to the deep water variance density spectra Sq^Oq) by 
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Figure 3.1 
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WIS Phase 1 Grid 
Locations Indicated. 
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MERCATOR PROJECTION OF PACIFIC OCEAN 



s h (f,0) = K 2 r (f,e) K 2 s (f) s 0 (f 0 ,e 0 ) J 



(3.4) 



where Kr is the refraction coefficient, Ks the shoaling coefficient, and J is the Jacobian 
of the transform. 

Here, J = <3(f o ,0 o )/5(f,0) 

however, since f Q = f, J becomes <3(0 O )/ 5(0) and integrates away. 

The shoaling coefficient, Ks, is given from Airy wave theory to first order as the 
ratio of the deep and shallow water group velocities. 

Ks = V(C g0 /C g ) (3.5) 

where C g is given by 

C g = C/2(l + 2kh/sinh 2kh) (3.6) 

In shallow water, surface gravity waves become non-dispersive and the group velocity 
can be approximated by 

C g = C = V(ghl . (3.7) 

The deep water group velocity is fixed for a particular frequency, whereas the shallow 
water group velocity will change with shoaling water depths; hence Ks will increase as 
the wave crest approaches the shore. 

The refraction coefficient, Kr, provides a measure of ray separation 

Kr = V(Vb) (3.8) 

where b Q is the deep water ray separation, and b is the shallow water ray separation. 
Values of Kr that are greater than unity represent energy convergence, and those less 
than unity, represent energy divergence. The refraction coefficients for this study were 
calculated using the linear refraction model by Dobson (1967) modified to refract rays 
from a specified shallow water location, back to deep water. This "back" refraction 
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model was developed to obtain refractive information at specific locations. The model 
searches for the 4-m isobath along a normal to the shoreline that originates at an 
observation site or other point of interest. Rays for each frequency in the WIS 
spectrum (0.03 Hz - 0.1 8Hz) are started at a specified location and depth, nominally 
4-m, and are propagated offshore at 0.1-degree increments over the range of all 
possible arrival angles. Once in deep water, the rays are stopped and the deep water 
angle is measured; then they are turned around and propagated back inshore along the 
same path to calculate the spectral refraction coefficient. The rays must return to the 
initial location within the specified area of ± 100 m to be a valid calculation. 
Examples of spectral refraction in Monterey Bay are presented in Figures 3.3, and 3.4. 

The bathymetry used in the refraction model is critical to the analysis. The 
accuracy of the calculated refraction information can be no better than the accuracy of 
the bathymetry; consequently, considerable effort was devoted to accurately depict the 
bottom. Original NOAA data was obtained that had been projected onto a six-second 
modified Universal Transverse Mercator (UTM) grid. The bathymetry data was 
initially screened for bad points, then the data was scanned along meridians and 
parallels for changes in slope that exceeded 30°. Points that generated unrealistic 
bathymetry were extracted from the data base, and the resultant bathymetry was 
projected onto an x-y plane via a modified UTM projection. The bathmetry was 
interpolated to a 200-m rectangular grid using a piecewise-linear, triangular plane 
interpolator (Franke, 1986). The triangulation method provided reasonable results in 
data sparce regions, and minimum distortion of bathymetric features. Intermediate 
smoothing was accomplished using a nine point weighted linear averager, and further 
smoothing was provided by the model itself which calculates bottom curvature by 
fitting a quadratic surface to adjacent isobaths. Waves originating from the northerly 
most quadrant are refracted when travelling over the shoal shelf region between the 
Farallon Islands, and Point Santa Cruz. To include the refraction to the north, a 
supplemental Northern bathymetry was gridded for a section of the California coast 
north of Santa Cruz. Refraction coefficients were calculated for the portion of the ray 
path which traversed the northern bathymetry, and for the portion of the ray path in 
Monterey Bay. The two refraction coefficients were multiplied together to determine a 
complete refraction coefficient for use in the spectral transformation. The Monterey 
Bay bathymetry is shown in Figure 3.5. 
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Figure 3.3 Ray Trace for Marina, f = 0.05 Hz.. 
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The objective in performing the spectral transformation is to produce a wave 
height that can be used to calculate wave set-up, run-up, and energy. The significant 
wave height, Hs, is related to the total energy, or variance, cr , by 

H $ = H rm$) ~ V2V(8"tf 2 ) ~ 4c (3.9) 

The variance is calculated by integrating the transformed shallow water spectrum over 
all frequencies and directions. 

= We s h (f, 0 )dfde (3.io) 

C. WAVE SET-UP AND RUN-UP 

As waves break in the surf zone, they transfer momentum shoreward. The 
change in shoreward momentum flux is primarily balanced by hydrostatic pressure; 
accordingly, the hydrostatic pressure balance is maintained through an increase in the 
mean sea level shoreward of the surf zone. This balance assumes that there are no 
mean currents, and no longshore gradients in the shoreward momentum flux, and is 
given by 

d/dx[ 2EC g /C - E/CJ = -pg(q + h)5q /dx ' (3.11) 

where the term in brackets on the left hand side is the momentum flux due to the 
normally incident breaking waves. As waves break, the wave energy, E, decreases 
resulting in a decrease in momentum flux, i.e. negative gradient, which must be 
balanced by a positively sloping sea surface. This mechanism is known as set-up (r|) 
and its maximum height above mean sea level occurs at its shoreward limit. Wave 
set-up studies by Guza and Thornton (1981) suggest that wave set-up is proportional 
to the deep water RMS wave height. Additional work by Holman and Sallenger (1985) 
over a wider range of wave and beach parameters, corroborate the results of Guza and 
Thornton, and further suggest that there is a correlation between the tidal stage, and 
wave set-up. Beach erosion primarily occurs at high tide, and consequently this study 
uses wave set-up at high tide (Holman and Sallenger, 1985). 
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Figure 3.4 Ray Trace for M.B. Academy, f = 0.09 Hz.. 
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(3.12) 



n - 0.88H' s0 

The significant wave height, H s is calculated at the 4-m isobath, where H $ ~ 4<r , then 
the wave is back shoaled without refraction to obtain a deep water significant wave 
height, H' $0 , using 

H' s0 = Hs/Ks (3.13) 

where Ks is calculated using the peak frequency. Wave set-up is dipicted in Figure 3.7. 

Wave run-up is the maximum superelevation of the sea level caused by the rush 
of breaking waves up the beach. Thornton and Guza (1981) suggest that run-up is also 
proportional to the deep water wave height H' s0 . Holman and Sallenger (1985) note a 
tidal correlation and find run-up at high tide to be given by 

n = 0.14H' S0 (3.14) 

Wave set-up and run-up are important parameters in raising the mean sea level to a 
stage where waves can erode the cliff base. 

D. TIDES 

The calculation of tidal height is inherently critical to the specification of mean 
sea level. The simultaneous presence of high tide and storm waves presents the highest 
probability that an erosion event will occur. Consequently, any parameterization of 
the fluxuating sea surface must include a harmonic tidal analysis. The harmonic 
analysis allows one to express the tidal height at any time, t, as a finite sum of its 
harmonic constituents, each of which is representative of a specific astronomic 
condition. The tidal expression from Dennis and Long (1971) is presented below 

h = H 0 + £ n f n H n cos[ a n t - (k n - [ V Q + u] n )] (3.15) 

where H Q is the height of the mean tidal level (above datum), f R is the node factor of 

the constituent, H n is the mean amplitude of the constituent, a n is the speed of the 

constituent, t is the time reckoned from some arbitrary point, k n is the epoch of the 

constituent at t = 0 for the observation period, [Vq + u] n is the value for the 
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argument at t = 0 for the observation period, and n is the particular constituent being 
computed. 

The constituents for Monterey Bay were obtained from a harmonic tidal analysis 
performed by the Tidal Predictions Branch of the National Ocean Service. Their data 
base was a thirteen-year time series of tidal data observed at Monterey's Municipal 
Wharf No. 2. Twenty constituents were isolated with the threshold criterion that a 
constituent would be consideried valid only if its amplitude contribution was greater 
than 0.6 cm. The speed, node factor, and argument of the constituents were extracted 
from Schureman (1971). The twenty constituents were used to calculate tidal heights 
at every hour during the twenty-year WIS time series. 

The tides in Monterey Bay are typically semi-diurnal with two low and two high 
water cycles per day. The mean tidal range is 0.98 m. 
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Figure 3.7 Wave Set-up. 
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IV. AN EROSION MODEL FOR MONTEREY BAY 



A. MODEL PRESENTATION 

The objective of this thesis is to develop an erosion model for Monterey Bay. 
Lima and Sklavidis (1984) computed a multiplicative constant, k, for the Brunn erosion 
model presented in Komar (1983) 

R = k[(S - T)/tanP] (4.1) 

where R is the recession rate, S is the height of mean sea level, T is the cliff toe 
elevation, tanP is the beach slope. The value k = 0.000098 was calculated from a 
linear regression of the Lima and Sklavidis (1984) recession rate data to the Brunn 
model. 

This research presents a model where recession is a non-linear function of the 
deep w'ater significant wave height H' s0 (for purposes of clarity in the following 
equations, H' s0 will be replaced by H). 

R mod = [AH 2 + BH + C(S - T)]/tanP (4.2) 

R mo d * s t ^ ie modeled recession rate and A, B, and C are constants. The deep water 
significant wave height is determined as discussed in Chapter III and the mean sea 
level, S, is given by 

S = Tide + Set-up + Run-up (4.3) 

The beach slope, tanP, was extracted from large scale hydrographic charts and varied 
little along the Bay. The idea of the model is that bluff recession only occurs when 
superelevation of the mean sea level, caused by tides and waves, exceeds the height of 
the cliff toe. When this criterion has been met, recession will occur at a rate which is 
proportional to: the wave energy, All 2 , the significant wave height, BII, and the 
ammount of w'ater which exceeds the cliff toe, C(S-T). 
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In solving for A, B, and C, it is necessary to minimize 



K R mod - R obs> 2 - MIN «-4> 

where R obs are the recession rates observed during the time-span of the WIS spectral 
data from 1956 to 1975. The partial derivatives of equation 4.4 are taken with respect 
to A, B, and C, and set equal to 0 to define the following three normal equations 

a£H 4 + BXH 3 + c£(S-T)H 2 = TanP lR obs H 2 (4.5) 



a£h 3 + b£H 2 + cT(S-T)H = Tanp £R Qbs H 



(4.6) 



a£H 2 (S-T) + B£H(S-T) + c£(S-T) 2 = TanP £R obs (S-T) (4.7) 



The three equations are then solved simultaneously by the matrix equation 



£H 4 £h 3 I(S-T)H 2 




’a" 




TanP lR obs H 2 


£h 3 £h 2 £(s-T)h 


X 


B 


zz 


TanP lR obs H 


£h 2 (s-T) ^h(s-t) I(s-t) 2 




C 




Tanp lR obs (S-T) 



This procedure is repeated for each of the twelve stations and the resultant coefficients 
are averaged to determine a mean A, B, and C for Monterey Bay. 

The model is very sensitive to changes in the elevation of the cliff toe, as this 
elevation is the model's go, no-go threshold value. An initial mean sea level model was 
run at each station to determine the number of hours, N, that the sea height exceeded 
the cliff toe for the twenty year WIS data. Adjustments , less than 0.5 m, were then 
made to the surveyed cliff toe elevation to ensure that the modeled N for each station 
was physically reasonable, and its value was within one order of magnitude of the other 
stations. The adjustments themselves are not unreasonable considering the inherent 
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variability of the cliff toe elevetion due to falling cliff material, and seasonal sediment 
transport. 

B. DISCUSSION OF RESULTS 

Recession model coefficients A, B, and C are determined for all sites that exhibit 
recession between 1956 and 1975. Moss Landing Marine Laboratory, and North Moss 
Landing are the only sites where continuous accretion during the time period prevented 
application of the model. Initial variance and standard deviation values are determined 
for the coefficients at the remaining ten sites, and based on this initial analysis, the 
coefficients for Monterey Sand Co. are rejected as being statistically unrepresentative 
of the erosion processes in Monterey Bay. The coefficients of the remaining nine sites 
are averaged to provide a mean A, B, and C, and hence define the recession model. 
These coefficients are listed in Table X. 

The calculated mean values of the coefficients provide a reasonable depiction of 
bluff recession for regions of low wave energy, and will heretofore be referred to as the 
Low Energy Model. The linear wave height coefficent B, is the most significant, and 
defines a mildly sloping recession profile. Unfortunately, as can be seen in Figure 4.1, 
the Low Energy Model becomes unreliable in regions of energy convergence. A 
significant wave height of 5 m appears to be the threshold where the model's energy 
term, AH 2 , begins to have a significant impact on the model recession. Hence a new 
value, Apjg , for the energy coefficient, A must be determined from a weighted average 
of the A coefficients computed in the following regions of high wave energy: South 
Fort Ord, Stillwell Hall, Marina, and Monterey Bay Academy (see Table X). A High 
Energy Model is then obtained by replacing A with Apjg . The coefficients B and C 
remain the same as in the Low Energy presentation, however, the new, Ajj£ 
parameterization of the modelled energy produces a much more realistic model 
response in energetic regions. This improved response is also depicted in Figure 4.1. 
The high energy recession curve is too shallow at low wave heights and consequently, 
the model overpredicts the observed recession in these regions. 
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TABLE X 

CALCULATED RECESSION COEFFICIENTS 


Station Name 


A [1/M 2 ] 


B [1/M] 


C [1/M] 


NPS Beach Lab 


-0.058619 


0.233661 


-1.79972 x 10’ 2 


Phillips Petroleum 


-0.025892 


0.144141 


-2.29962 x 10‘ 3 


Sand City 


-0.052943 


0.358974 


-2.59265 x 10‘ 3 


Monterey Sand Co. 


-0.063531 


0.461977 


-1.83579 x 10‘ 2 


South Fort Ord * 


-0.011144 


0.133575 


5.18795 x 10' 4 


Stillwell Hall * 


-0.023181 


0.233329 


2.43953 x 10’ 3 


Marina * 


-0.016920 


0.196066 


-1.20543 x 10' 3 


Rincon Beach 


-0.010737 


0.101922 


5.60665 x 10* 5 


Zmudowski Dunes 


-0.020776 


0.188279 


-2.08930 x 10’ 3 


Monterey Academy * 


-0.023647 


0.228347 


3.02067 x 10’ 3 


Mean Values 


-0.027095 


0.20203 


-0.002317 


High Energy values * 


-0.020051 


same 


same 


Standard Deviations ±0.0287 ±0.2143 

* Indicates regions of high wave energy, H' s0 ^ 5 m. 


±0.0024 
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Figure 4.1 Recession Model Profile. 
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RECESSION (METERS) 



C. CONCLUSIONS 



The recession rates are found to be related to the local wave height. A 
polynomial regression model in wave height shows the recession to be most correlated 
with wave height to the first power (B is the largest coefficient) and negatively 
correlated with wave energy, H 2 . The recession observations at Sand City and 
Monterey Sand Company show significantly more bluff recession than the Low Energy 
Model suggests, and the observed significant wave height does not warrent use of the 
High Energy Model. Local sand mining operations are the most probable cause for 
the high recession values at these sites. For the remaining eight stations, however, the 
most representative results are obtained by applying the Low Energy Model in regions 
of low wave energy, and the High Energy Model in regions of high wave energy. In 
doing so, the standard error calculated between the observed and the computed 
twenty-year recession rates is ±3.41 m which equates to ±0.17 m/yr. The mean 
recession rate for the stations which experienced recession between 1956 and 1975 , 
discounting stations Sand City and Monterey Sand Company, is 0.76 m/yr, hence the 
model represents actual recession to ±22 %. These results are listed in Table XI. 
Graphic depictions of the twenty year recession at each site are depicted in Appendix 
C. The observed recession results presented in Table XI only represent events occuring 
between January 1956 and December 1975, and a complete summary of bluff recession 
between 1940 and 1984 is provided in Table XII. 



70 



TABLE XI 

RECESSION MODEL RESULTS (1956-1975) 



Station Name Observed Modelled High Energy Model 



NPS Beach Lab 


9.00 m 


11.45 m 




Phillips Petroleum 


7.80 m 


13.43 m 




Sand City 


23.80 m 


14.4 m 




Monterey Sand Co. 


33.0 m 


14.86 m 




South Fort Ord 


15.8 m 


10.57 m 


19.68 m 


Stillwell Hall 


19.40 m 


10.97 m 


20.49 m 


Marina 


22.4 m 


12.36 m 


23.49 m 


Rincon Beach 


9.60 m 


13.56 m 




Zmudowski Dunes 


16.80 m 


14.08 m 




Monterey Bay Academy 


21.8 m 


13.51 m 


19.97 m 


The Standard Deviation for the Model is ± 
years which equates to ± 0.17 m/yr. 


3.41 m for 20 
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TABLE XII 

OBSERVED RECESSION SUMMARY 



Station Name 


Recession (m/yr) 


NPS Beach Lab 


-0.58 


Phillips Petroleum 


-0.85 


Sand City 


- 0.97 


Monterey Sand Co. 


-1.94 


South Fort Ord 


-1.12 


Stillwell Hall 


-1.98 


Marina 


i 

O 
L j 
oo 


Rincon Beach 


-0.84 


Marine Lab 


+ 0.45 


North Moss Landing 


+ 0.03 


Zmudowski Dunes 


-0.52 


Monterey Bay Academy 


-0.33 


Standard Deviation 


±0.73 Mean = -0.75 m/yr 
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APPENDIX A 

PROGRAM LISTING PHOTODAT 



— I N*- , U' r "WHAT IS VOUR ST mT I QN I. D. I. 5. LAB1336" ;E$ 

3 OPEN £3 FOR OUTPUT PS *1 

4 -'FIN 7 

5 PRINT NILE'S 

7 PRINT 

8 CG(JNT=i3 

3 PRINT " REMEMBER “0 ALWAYS APPROACH TOUR POINT FROM THE SAME DIRECTION 1 ' 

10 PRINT "TO PREVENT ELECTRONIC OR MECHAN ICPL BACKLASH r ROM OCCURING IN" 

11 PRINT "THE DATA. THE -IRST STEP IS “0 DEFINE YOUR ORIGIN. (MOVE THE" 
13 PRINT "VIEwER OVER THE CORRECT POINT PND PRESS ( DH ) (RETURN) PT THE" 

13 PRINT " ( : ) PROMPT. RE D EPT “HIS PROCEDURE FOR EPCH AXIS." 

14 PRINT 



15 GOTO 30 
IE PRINT 
18 PRINT 

13 POINT " REDE- I NE THE ORIGIN, PRESS ( 
30 FOR 1=1 TO 3 STEP 3 
33 J=I -1 
30 \=1000+J 
33 P= : 

40 PRINT CHR$ < INP <N) ) ; 

4? B=I NP < N+l ) 

50 C= I NT < 5/3 ) 

33 IF C=B/3 THEN 30 

£0 IF I NT (C/3 ) =C/3 THEN 45 

£5 IF A=4 THEN 33 

70 AS=INKEY* 

71 IF i_EN(P3)=0 THEN 45 
73 PRINT PS; 

73 CUT ASC(PS) 

77 P=P -1 

78 GOT 0 45 

80 P$=CHRs ( INP (N ; .< 

83 IF’ AS <> CHRS ( 10) THEN PRINT PS; 

30 IF Ps=CHRS( 63> “HEN 18 



Dn ) PT THE PROMPT" 

’ THIS LOGP INITIALIZES THE XSY AXES 
’ N= THE DMC-400 READ /'HR I TE ADDRESS 
’ N = 1000 X-AXIS N= 1003 Y-PXIS 
’ AXIS COUNTER. X = 1 Y=3 

’ CHECKS DMC-400 STATUS ADDRESS 
’ "C" CHECKS STATUS FLAG BIT 

* DOES DMC-400 HAVE DATA TO SEND n 
' IS DMC-400 READY TO ACCEPT . COMMANDS 

’ INPUT COMMAND FROM KEYBOARD 



1 SEND THE "DEFINE HOME" COMMAND 



’ READ DATA FROM THE DMC-400 



34 GOTO 43 

35 NEX“ I 

100 PRINT "LEAVE Tr-E X/Y CONTROLS ALONE AND CHECK TO INSURE TH A T X & Y =0 " 

103 PR I N“ 

105 INPUT "PRESS 1 p ) WHEN READY" ;FS 
110 IF r $=' P " .THEN 133 
115 GOTO 105 
133 PRINT 

130 INPUT "HAS “*-£ ORIGIN BEEN CORRECTLY DEFINED"' (Y/N)";DS 
133 IF Ds= " Y " THEN 143 
140 GOTO 1 £ 

14<* PRINT 

145 PRIN“ "YOU APE NCw READY TO BEGIN TAKING DATA. MOVE THE VIEWER OVER" 

14B PRINT "THE DESIRED POSITION AND PRESS ( DH ) TO CATALOG T h-E POINT. START 

147 PRINT "WIT— ~V£ BASELINE, THEN SCALING POINTS. THEN THE SHORELINE POINTS. 

148 PR IN- 
133 1 = 1 

300 “OR 1=1 TO 5 STEP 3 ’ LOOP TO TAKE X/Y DATA 

310 J=I-1 
330 N=:000+j 
333 PRINT 

335 PPIN- "-'RESE ' “P ) ”0 C A — A ■_ C 3 " ; N : " ADDRESS DA T A. " 

330 GCSLB 500 ’ SUBROUTINE TO CO 

50 I “ N=: 00 i? ~-E' 33/ 



MMUNICATE W/DMC -^00 



26 0 X=?-l)*9 

223 5CTQ 278 

177 IF N - 1 002 ~HEN 171 

171 GCTO £75 

272 Y=R*(-1) 

272 GOTO 278 

275 Z=R 

276 NEXT I 

273 ='RINT 

26*21 PRINT "YOUR X, Y & Z COORD’S ARE ";X,Y, Z 

281 INPUT "WOULD YOU LIKE TO STORE THESE PTS IN THE DATA FILE (Y/N):":YS 

282 IF YS=*"N" THEN 283 

283 PRINT 41, X ;Y;Z 

284 COUNT =COUNT 1 

285 IF X=0 THEN 123 

286 IF Y=0 THEN 123 
283 PRINT 

230 INPUT "WOULD YOU LIKE TO CATALOG MORE DATA POINTS - ’ (Y/N)";G$ 

300 IF GSO'VN" THEN 148 

302 PRINT 

203 PRINT "THERE ARE COUNT ; " POINTS STORED IN FILE ";ES 
304 PRINT 

205 PRINT "BE SURE TO LOG THIS INFORMATION, AS YOU WILL NEED IT LATER 1 " 
306 PRINT 

308 PRINT "THAT’S ALL FOLKS!!' ITS MILLER TIME !!!" 

303 CLOSE 
310 END 

500 REM THIS IS A SUBROUTINE TO COMMUNICATE BETWEEN THE DMC-400 S IBM PC 
502 A= 1 
505 M=6 
510 R=0 
515 K= 1 

520 PRINT CHRS ( INP (N) ) ; 

530 B- I NP ( N-*- 1 ) 

540 C- I NT ( B /2 ) 

550 IF C=B/2 THEN 565 

555 IF INT (C/2) =C/2 THEN 520 

556 IF A=4 THEN 630 

557 AS = INKEYS 

558 IF LEN ( AS ) =0 THEN 530 
553 PRINT AS; 

560 OUT N, ASC(AS) 

562 A=A+ 1 

564 GOTO 530 

565 AS=CHRS ( INP (N) ) 

566 IF AS <) CHRS (10) THEN PRINT AS; 

568 IF AS=CHRS (63) THEN 502 
570 V=ASC (AS) 

580 IF (V>=48) AND ( V < = 57) THEN 620 
530 IF (V>=65) AND ( V < = 70) THEN 620 
610 GOTO 530 

620 IF V (60 THEN Z=(V-48) 

630 IF V) 60 THEN Z= < V-55) 

64 0 R=R + ( Z* 16 "- (M-l > ) 

645 IF M=6 THEN Z2=Z 
650 M=M- 1 

660 IF K=6 THEN 630 
£70 K=K+1 
680 GCTO 530 
630 IF Z2 < 8 THEN 710 
700 R=R— 1 67772 164 
710 RETURN 



\ THIS SUBROUTINE COMMANDS THE DMC-400 
’ TO REPORT ITS POSITION (TELL POSITION 
’ COMMAND), AND THEN READS THE POSITION 
’ FROM THE N ADDRESS- THE POSITION (IN 



’ CHECKING STATUS & READING DATA 
’ AXIS COUNTER 1=X 2=Y 

’ ASSIGN "TELL POSITION COMMAND" 
’ SEND COMMAND TO DMC-400 



’ CHECKING TO SEE IF DATA ARE 
’ CHARACTERS 0-3 OR A-F 

’ CHANGING FROM HEX TO DEC 
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APPENDIX B 



PROGRAM LISTING WRKPROG 



10 INPUT "INPUT THE NUMBER OF POINTS IN THE DATA FILE:":N 
£0 DIM X# (N) . Y# (N) , Z# (N) , XR# (N) , YR# (N) 

30 PRINT 

40 INPUT "INPUT THE NAME OF YOUR DATA FILE, I.E. LAB19 33 : " ; E$ 
45 PRINT 

50 OPEN "I", #1,ES 

70 FOR 1=1 TO N 

00 INPUT #1, X# (I) , Y# (I) , Z# (I) 

90 NEXT I 
90 1=1 



100 AX# = . 000123762075# 

105 BX#=-. 00072479240# 

110 qy* = . 000093559173** 

115 BY#=-. 00513229147# 

150 FOR 1=2 TO N 
100 X# ( I ) = X# ( I ) *AX# + BX# 
170 Y# ( I ) = Y# ( I ) *AY# + BY# 
190 NEXT I 



REGRESSION COEFFICIENTS TO CONVERT 
ENCODER COUNTS TO CENTIMETERS. WHERE 
CM = A* ( E— COUNTS ) + B 



’ CONVERSION E— COUNTS TO CM 



200 ALFA#=ATN < Y# (2) /X# (2) ) * COORDINATE SYS ROTATION ANGLE 

210 FOR J=1 TO N 

220 XR# (J) « (X# (J) *COS (ALFA#) ) + (Y# (J) *SIN (ALFA#) ) ’ COORDINATE TRANSFORM AT I ON 

230 YR# (J) = ( (-1 ) *X# (J) *SIN (ALFA#) ) + ( Y# ( J ) *COS (ALFA#) ) 

240 NEXT J 
245 PRINT 

248 INPUT "INPUT THE MEASURED FIELD DISTANCE (IN METERS ) : " ; GD I ST 
250 PRINT 

270 PDIST#=SGR ( (XR# (4) -XR# (3) > ^2 + (YR# (4 ) -YR# (3) )^2) 

274 LPRINT 



275 LPRINT "PHOTOD I STANCE = " ;PDIST#;" CM 
277 PRINT 

200 PSCALE#= (PD I ST#/ 100) /GDI ST 
233 PRINT PSCALE# 

205 P«=i /PSCALE# 



239 LPRINT 



290 LPRINT "PHOTOSCALE IS: " ; P# 

291 LPRINT 

292 LPRINT "YOUR ACTUAL SCALE, ROTATED X/Y PAIRS FOLLOW (METERS) 
300 FOR K= 5 TO N 

310 X R # ( K ) = (XR# (K) /PSCALE#) / 1 00 
320 YR# (K ) = (YR# (K) /PSCALE#) / 100 
335 LPRINT XR# ( K ) , YR# (K) 

350 NEXT K 
360 LPRINT 



400 FOR L = 5 TO ( N— 1 ) ’ TRAPAZO I DAL RULE INTEGRATION 

420 AREA#=AREA# + ( (YR#(L)+YR#(L+1) )/2) *ABS (XR#(L+1)-XR#(L) ) 

450 NEXT L 
460 PRINT 

500 LPRINT "THE AREA IS: " ; AREA# ; " SQUARE METERS FOR STATION " ; E$ 

600 CLOSE 
900 END 
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20YR RECESSION (METERS) 
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20YR RECESSION (METERS) 
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20YR RECESSION (NETERS) 
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